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We have determined families of two-dimensional deterministic totalistic cellular automa- 
ton rules whose stationary density of active sites exhibits a period two in time. Each 
family of deterministic rules is characterized by an "average probabilistic totalistic rule" 
exhibiting the same periodic behavior. 

Many natural populations of plants and animals exhibit large fluctuations of 
density with a roughly cyclic behavior. A well-known example is the oscillatory 
behavior of the Canadian lynx population as documented in the data compiled by 
the Hudson Bay Company over the period 1735-1940. Oscillations with an approxi- 
mate period of about 10 years are observed with lara* amplitude fluctuations, which 
could, actually, correspond to a chaotift behavior.liH3 jGood introductions to popu- 
lation dynamics may be found in MayEI and Murray.Q Most models in population 
dynamics are formulated in terms of differential equations or difference equations, 
which means that the local character of the interactions between prey and preda- 
tors, for example, is not taken into account. In order to describe more correctly the 
local character of the predation process, it would be better to formulate predator- 
prey models in terms of cellular automata (CAs), which are fully discrete dynamical 
systems evolving in time according to local rules. This paper does not deal with 
any specific model. Its purpose is to understand undoj which circumstances CAs 
exhibit nontrivial global behaviors. It has been arguedB that, in spatially isotropic 
systems with short-range interactions evolving in discrete time, periodic states with 
periods larger than 2, are never ratable under generic conditions. Despite these 
arguments, Chate and MannevilleEHa gave several examples of high-dimensional de- 
terministic CAs exhibiting a global (quasi)period-3 behavior^. In order to verify that 
the corresponding states were not metastable Gallas et alu performed a detailed 
study of the Chate-Manneville CAs, and confirmed that the nontrivial behavior 
found by these authors extends to much longer times and to much larger lattices. 
Since the chance to find such a behavior should increase with the dimension of the 
CA, it is of interest to find low-dimensional CAs exhibiting the same behavior. A 
three-dimensional deterministic CA exhibiting a quasiperiodic behavior has been 
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found by Hcmmingsont3, and shortly after its nontrivial cpUective behavior has 
been confirmed through large-scale simulations (Chate et aLtil). 

In this note we present families of two-dimensional totalistic class-3 CAs which 
exhibit a period-2 behavior. The class-3 requirement eliminates CAs whose attrac- 
tor consists of only two points in the configuration space. No arguments have been 
given against such a collective behavior, but these systems might be of interest since 
they are not so frequent. Moreover, for each neighborhood defining a CA family, 
there exists a well-defined "average probabilistic CA rule" which exhibits the same 
collective behavior. 

Two-dimensional deterministic totalistic CAs may be defined as follows: Let 
s:ZxZxNi— >{0, 1}bea function that satisfies the equation 

s(zi, j2,t + 1) = /( ^ s{ii+ ji,i2+ j2,t)j, 
{jij2}ev, 

for all (ii,?2) G and all i e N, where Z is the set of all integers and N the set 
of nonnegative integers. Vr defines a neighborhood of site {11,12) which is either of 
the "von Neuman-type" with 

Iii| + |j2| <r 

or of the "Moore-type" with 

|ji| < r and \j2\ < r 

The integer r characterizes the range of the neighborhood. For a given neighbor- 
hood, each totalistic rule is, therefore, determined by a function 

S>^f{S), 

where 

5g {0,I,2,...,5,nax} and /(5)e{0,l}, 

with S'niax ~ (2r + f )^ for a Moore-type neighborhood and S'max = 2r^ + 2r + 1 for 
a von Neuman-type neighborhood. 

For a given neighborhood, the number of rules is 2"^"-'"'. For relatively small 
neighborhoods S'max < 13, a systematic search for rules exhibiting nontrivial collec- 
tive behavior is possible. For the simplest Moore and von Neuman neighborhoods, 
with r = 1 we did not find any nontrivial behavior. The smallest neighborhood 
for which such a behavior has been found is the von Neuman neighborhood with 
r = 2, and the corresponding family of 2^"^ = 16384 rules has been systematically 
investigated. Systematic investigations of families of CA rules with range larger 
than 2 is impossible. In this case, for each family, we studied a randomly selected 
sample of about 50000 rules. 

von Neuman neighborhood with r = 2. 

Each rule has been first studied on a 128 x 128-lattice over 200 time steps. The 
average concentrations Cg and Cq of nonzero sites over the 50 last even and odd time 
steps were evaluated. We only retained rules satisfying the condition |ce — Co| > 0.1. 
A few hundred rules satisfying this criterion have been thoroughly investigated on 
a 512 X 512-lattice over 20000 time steps . Several behaviors have been observed. 

• Some of them are "similar" to Conway's Game of life, in the sense that only 
periodic small structures, with various periods, remain on a period-2 back- 
ground (see Figure 1). In this case, the attractor consists of a very few number 
of points in the configuration space. 



Periodic behavior in 2D Cellular Automata 3 




Fig. 1. Rule 189: periodic structures on a period-2 background. 



• Some converge to a fixed number of nonzero sites after large transient time 
during which the concentration oscillates. Their patterns are inhomogeneous 
(see Figure 2). The "local" concentration of nonzero sites oscillates; however, 
at large times, in the infinite-size limit, the concentration of nonzero sites tends 
to a fixed point (see Figure 2). Suchj-a. behavior has already been observed in 
some probabilistic cellular automatall3. 

• The 54 rules listed in Table I have a nontrivial collective period-2 behavior 
for the concentration of nonzero sites. They appear by pairs, rule g defined 

by 

giS) = 1 - f{S,na. - S) 

has the same properties as /. A typical example (Rule 2077) is shown in 
Figures 4 and 5. Note that all these rules are such that /(O) = 1. 

In order to characterize a family of rules {fk \ k e K}^ which exhibit a period 
two in the infinite-time limit, we define the function /average ■ [0, 1] >— > [0, 1] by 

/average ('S'/ Smax] TjFl ^ ^ fki^S)^ 

where \K\ denotes the number of rules belonging to the K family (see Figure 6). 

Rules with larger neighborhoods. 

We have studied three other families characterized by the following neighbor- 
hoods. 

• Range-3 von Neuman neighborhood: Sj^a^ — 25, i.e. 2^^ — 67108864 rules. 

• Range-2 Moore neighborhood: Smax = 25, i.e. 2^6 = 67108864 rules. 

• Range-4 Moore neighborhood: S'^ax = 81, i.e. 2^^ w 4.836 10'^^ rules. 
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Fig. 2. Rule 2158: evolution of the concentration c of nonzero sites on a 512 X 512-lattice. 





Fig. 3. Spatial patterns of Rule 2158 at two consecutive time steps after 2 x 10* iterations. The 
lattice has 512 X 512 sites. There are short range correlations over a length scale which is much 
larger than the range r = 2 of the rule. However the concentration of nonzero sites tends to a 
fixed point on an infinite lattice. 
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Fig. 4. Rule 2077: evolution of the concentration of nonzero sites as a function of time on a 
512 X 512-lattice. A period-2 behavior is clearly observed. 



Fig. 5. Patterns of the Rule 2077 at two consecutive times after 2 x 10^ iterations. 
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Tabic 1. Deterministic totalistic CA rules exhibiting period-2 collective behavior with a range-2 
von-Neuman neighborhood. 
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Fig. 6. Graph of the function /average characterizing the family of rules defined by a range-2 von 
Neuman neighborhood. 
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As mentioned above, a systematic investigation of such a large number of rules 
is impossible. Therefore, for each family, we only studied a sample of about 50000 
rules selected at random. We found that a small proportion of these samples did 
exhibit a period-2 behavior, i.e. about 0.7% of the rules with a neighborhood 
containing 25 sites, and 1% of the rules with a neighborhood containing 81 sites. In 
contrast with the family of rules with a range-2 von Neuman neighborhood, about 
40 to 50% of them were such that /(O) = 0. For these three families, we determined 
the functions /average whose graphs are represented in Figure 7. These functions 
are well approximated by a simple combination of two sine functions of the form: 

/averagc(a:) = 0.5 - ^sin(27r(x - 0.5)) - Bsin(27r(a; - 0.5)) 

where x — S/Smax- Fitting the data, wc foimd: 

• for the range-3 von Neuman neighborhood, A « 0.291, B « 0.008; 

• for the range-2 Moore neighborhood, A « 0.304, B « 0.003; 

• for the range-4 Moore neighborhood, A « 0.138, B « 0.074. 
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Fig. 7. Graph of functions /average characterizing famihcs of rules defined by: (a) range-3 von 
Neuman neighborhood (diamonds and dashed hne), (b) rangc-2 Moore neighborhood (open circles 
and solid line), and (c) range-4 Moore neighborhood (bullets and solid line). 



To a function /average we associate a probabilistic totalistic law defined by 

f(S) = I ^' probability p = faverage{S/ Sma;K); 

•'^ ' \ 0, with probability 1 - p. 
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Probabilistic totalistic CAs evolving in time according to rules of this type exhibit 
period-2 behaviors. As an example, the evolution of the concentration c of nonzero 
sites for the probabilistic rule associated to the family defined by the range-4 Moore 
neighborhood is represented in Figure 8. The concentration oscillates between ci = 
0.337 and C2 = 0.663. The mean-field approximation predicts the amplitude of the 
oscillations with a very good accuracy. Within this approximation, it is found that 
the concentration oscillates between c'^^'^ = 0.335 and c™^** = 0.665. 
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Fig. 8. Concentration of nonzero sites of a probabilistic totalistic CA whose rule is determined by 
the function /average characterizing the family of rules defined by the range-4 Moore neighborhood. 

The proportion of deterministic totalistic class-3 CAs whose density of nonzero 
sites is periodic in time increases with the range of the neighborhood characterizing 
the evolution rule. In agreement with an argument due to Grinstein,El no period 
larger than 2 has been observed. To each family of rules exhibiting a period-2 
behavior, it is possible to associate a well-defined mapping from [0, 1] to [0, 1]. Such 
a mapping defines a probabilistic totalistic CA rule, and the density of nonzero sites 
of a CA evolving according to such a rule also exhibits a period-2 behavior. For 
these probabilistic CAs, the mean-field approximation predicts the values of the 
density with a good accuracy. 
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